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ABSTRACT 

We adopt a model independent method to reconstruct the dark energy equation of state by analyzing 5 sets 
of SNe la data along with Baryon Acoustic Oscillation (BAO) and Observational Hubble Data (OHD). The 
SNe la data sets include the most recent UNION2 data and other data compilations from the year 2007 to the 
present. We assume a closed form parametrization of the luminosity distance in terms of redshift and perform a 
analysis of the observational data. The matter density at the present epoch iljjj is also taken to be a parameter 
in the analysis and its best-fit values are obtained for each of the data sets. We found a strong dependence of 
dark energy equation of state on the matter density in the present and earlier epoch. From the analysis, we also 
predict the lower limit of matter density parameter at an earlier epoch within Icr confidence level for a flat FRW 
universe. The dark energy equation of state appears to be a slow varying function of z- The variation of dark 
energy density parameter and the matter density parameter are also shown along with their Icr variations. 
Subject headings: Dark energy. Supernova la, Baryon Acoustic Oscillations, Observational Hubble Data 



1. INTRODUCTION 

Observations on type la supemo vae (SNe la) 
dRiessef g/. 19981: iPerlmutter ef a/. 1999h during nearly 
last two decades reveal that the universe is undergoing 
accelerated expansion in the present epoch. This accelerated 
expansion can be explained by invoking the existence of dark 
energy - a hypothetical energy component with a negative 
pressure. Despite several past and ongoing efforts the nature 
and origin of 'dark energy' remains a mystery. The evidence 
of dark energy can also be predicted from th e observation of 
Baryon Acoustic Oscillation (BAO) (lEisenste in et al. 2005]), 
Hubble data based on differential ag es of the galaxies (OHD) 
(iJimenez et al. 20021: iMa et al. 201 II) etc. 

Some of the properties of dark energy can be extracted by 
performing analysis of the observational data. In general, 
there are different approaches for the analysis of SNe la 
data (obtained in the form of luminosity distance modulus 
versus redshift). One of the approaches involves choice of 
some arbitrary parametrization of the dark energy equation 
of state wx(z) (=p/p). p and p respectively denoting the 
pressure and energy density associated with the dark energy. 
The luminosity distance di(z) and the bolometric magnitude 
/i(z) at redshift z are then found by using the assumed 
parametric form of wx{z). Such approaches are taken 



up and discussed in detail in 


Starobinskv ef flZ. 1998 


Huterer et al. \ 9991 lAstier 2000 


: Chiba et al. 2000 


WeUer et al. 2002 


: iMaor et al. 2002 


IChevallier ef fl/. 2001 


Linder et al. 2003 


; Jassal et al. 20051 


Sahni et al. 2003 


Gerke et al. 2002 


: ICorasaniti et al. 20031 


Wang et al. 2004 
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i Nesseris et al. 2004t I Roy Choudhurv et al. 2005 : 

ICong 2005^, JWetterich 2004 IWuef a/. 2006 : 

Guo et al. 20051: ISimon et al. 2005?). In another kind of 
approach jiiz) are first fitted with observational data and 
then one finds the dark ene rgy e quat ion of state wxiz) 
(iDalv et al 200l F av et al . 2006t iHuterer et al. 2005t 
ISaini 2003H Shafieloo 2007c i Wang et al. 20011) . 

In this paper we have considered differ- 
ent compilations of SNe la data sets viz. 

fess et al. 2Q01\ I WoodVasev ef al. 20071: iDavis et al. 20071: 
A'alski et al. 2008' iKessIer et al. 20091) aiid 

nanullah et al. 20Tol) . From SNe la observations these 
groups have tabulated the values of /i(z) for different values 
of redshift (z) within the Hmit 0.001 < z < 1.76. We have 
considered a parametric form of diiz) and expressed jiiz) in 
terms of this parametric form of diiz)- We make a ^ analysis 
of the combined data sets of SNe la, BAO and OHD to obtain 
the best-fit values of the parameters from the observational 
data. We have taken 5 different SNe la data sets (described 
later) and for the combined analysis each of these data sets 
are combined with BAO and OHD data. In this work the 
matter density at the present epoch (^1°) ^1^° taken to be 
a parameter for the \^ analysis, and by performing the 
minimization we obtain the best-fit value of (fij],) along with 
the parameters appearing in the parametrization of diiz). 
Having thus obtained d^iz) as a function of z we calculate the 
variation of LOxiz) as a function of z for the best-fit values of 
the parameters and their Icr limits as well for each of the 5 
data sets considered. In all the calculations we assume that 
present universe is spatially flat and contains only matter and 
dark energy. The results of the analysis show that knowledge 
of the matter density of the universe at some earlier epoch is 
instrumental in providing observational evidences in favour 
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of varying dark energy or cosmological constant solutions. 
We have also shown the simultaneous variation of matter 
density parameter rim(z) and dark energy density parameter 
flxiz) with z for the best-fit values of the parameters (obtained 
from fitting) and their Icr range. We also found the epoch 
at which the dark energy started dominating over the matter 
component of the universe. 

The paper is organized as follows. In Section |2] the for- 
malism for reconstruction of dark energy is described. The 
methodology of the analysis are given in Section |3] In Sec- 
tion|4]we describe results of analysis of different sets of data 
and explain our results. Finally in Section |5] we make some 
concluding remarks. 

2. RECONSTRUCTION OF THE EQUATION OF STATE OF DARK 
ENERGY 

In Standard FRW cosmology, for a spatially flat universe, 
the luminosity distance di{z) of an object at a redshift z is 
related to the Hubble parameter H{z) as 



(1) 



where c is the velocity of light and a is the scale factor, whose 
time evolution determines the time evolution of the universe. 
Modeling the total content (other than gravitational field) of 
the universe as a perfect fluid characterized by its energy den- 
sity p and pressure p, the (00) component of the Einstein's 
equation for a spatially flat universe (K = 0) gives 
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Observations from WMAP experiment suggest that the 
present universe is spatially flat and contribution of ra- 
diation to the total density of the universe in negligible 
jKomatsu et al. 20Tll) . Measured values of redshifts of Su- 
pernova la events (z < 2) correspond to epochs close to the 
present epoch in cosmological time scale. Therefore for anal- 
ysis of Supernova la observations we are permitted to use Eqs . 
([T]i and dU. In this context we can also write the total en- 
ergy density of the universe by neglecting the radiation energy 
density as p{z) = Pmiz) + px(z) where p,„(z) and px(z) denote 
the contribution in the energy density due to (non-relativistic) 
matter and dark energy. This can also be expressed in terms 
of the corresponding density parameters H.^ = (3i/^/87rG)p,„ 
and ilx = OH^/^T^G)px as r2,„(z) + r2x(z) = 1- Assuming 
an effective equation of state for the dark energy wxiz) = 
Pxiz)/ px(z), the Hubble parameter can be expressed as 



H\z) 



= nO(l+z)3 + fl°exp( / 3(1 + wxiz!)) 



dz' 
l+z' 



,(3) 



where //q is the value of the Hubble parameter at the present 
epoch. and fix respectively denote values of matter and 
dark energy density parameter at the present epoch. Differ- 
entiating both sides of equation Eq. (|3]l with respect to z we 
obtain the equation of state of dark energy as 



wx(z) = -l + 



2 (l+z) dH(z) qO 

3 W(z) dz """H^feVHj 



(l+z)' 
"'HHz)/H' 



(4) 



We use the above equation to reconstruct the dark energy 
equation of state from observational data. Such a recon- 
struction therefore requires extraction of the quantities H(z) 



and from the observational data. From the measured 
values of the luminosity distances {di{z)) of type la Super- 
novae at different redshifts (z) we can obtain H{z) at dif- 
ferent redshifts using Eq. ([T]i. The observational Hubble 
Data (OHD) based on the differential ages of the galaxies 
also provide values of the Hubble parameter H(z) at some 
redshift z values. On the other hand, the measurement of 
Baryon Acoustic Oscillations (BAO) from the study of large- 
scale correlation function of sky surveys of several thousands 
of luminous red galaxies provides a val ue for the quantity 
jEisenstein et al. 20051; IZhang et al. 20101) 
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A(zi)=- 



[H{zi)/Ho\ 



1/3 



1 



dz 



ZiJo Hiz)/Ho 



(5) 



with zi = 0.35. Eq. Q thus connects the H{z) and 0° to the 
observed quantity A(z = 0.35). A joint analysis of the SNe la, 
OHD and BAO data can therefore simultaneously constrain 
the parameters H(z) and fljjj. 

The Hubble parameter H{z) and its variation with z show up 
in different ways in the SNe la, OHD and BAO observations. 
The Hubble parameter is again related to the luminosity dis- 
tance, diiz) by Eq. ([T]i. Choice of some parametric form of 
diiz) therefore relates the observed quantities with the param- 
eters of diiz). In this work we use the parametrization of the 
luminosity distance as 



dLia,b;z) = 



Ho 



z(l+az) 



l+bz 



(6) 



which respect the conditions that di = at z = and dio^z for 
large z corresponding to the radiation dominated era. SNe la 
data are available for z ^ 1 .76 - an era which is dominated by 
matter and dark energy. Using this parametric form and Eq. 
([T]| we can express the Hubble parameter and its z-derivative 
in terms of parameters a and b as 

{\+zf{\+bzf 



H{a,b;z) = Ho 



1 + 2az+ {a — b+ab)z 



(7) 



The expression for H{a,b;z) thus obtained can be used to ex- 
press the quantity A (in Eq. (|5]l) in terms of a and b as 



A(a,b,nl)=- 



[H(a,b;zi)/Ho] 



1/3 



dz 



H(a,b;z)/HQ 



2/3 

(8) 



with zi = 0.35, where ilj],, along with a and b, will be treated 
as parameters to be determined from observation. The method 
of analysis of observational data for obtaining the parameters 
a, b and Oq is discussed in Sec. [3] Finally the dark energy 
equation of state w(z) as given in Eq. (|4]l can be expressed in 
terms of the parameters a, b and J^jJ, as 



wx(a,b,nl-z) = 



3(l+fc) 



l+b + 2bz- 



(l+z)(l+fcz)(n+(gfc+a-fc)z) 
l+2az+(ab+a-b)z^ 



1 



QO [l+2az+(flfr+g-fc)z^]^ 
"1 (l+z)(l+fe)'' 



(9) 

The best-fit values of the parameters a, b and fij), obtained 
from the analysis of the experimental data can be used in Eq. 
^ to obtain the variation of the equation state of dark energy 
with redshift. 

Also we note that Eq. (O can be identified with the equation 

n„,(z) + flx(z)=l with 

Hin^(l+z)' 



^miz)-- 



HHz) 



(10) 
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so that ^lx(z) = 



H\z) 



exp(/„'3(l+wx(z))f|) = l- 



^m{z)- The matter density parameter ri,„(z) as given in Eq. 
(fTOb can be written in terms of the parameters a, b and $7^, 
using the expression for H{a,b;z) as 



nr A oo ^ oo (1+^) (l+2fl^+(aZ7 + fl- 



2\2 



[(l+z)(l+fe)]4 ^^^^ 

With best-fit values of the parameters a, b and obtained 
from analysis of observational data we use Eq. (fllT i to com- 
pute the z-variation of the density parameters both for matter 
{VLm) and dark energy (^lx = 1 - f^,,,). 

3. METHODOLOGY OF DATA ANALYSIS 

The SNe la data remains the key observational ingredient 
in determining cosmological parameters related to dark en- 
ergy. In this paper we have considered different compila- 
tions of SNe la observation s including the recent UNION2 
data (lAmanullah ef g/. 2010I). T he other SNe la data sets con- 
sidered here are ( Riess et al. 2 001: WoodVasev et al. 2007; 
iDavis et al. 20071: iKowalski et al. 20081: iKessler et al. 2009). 



These different groups tabulated the values of the distance 
modulus fiiz) for different values of the redshift z from the 
SNe la observations. The distance modulus /i is related to the 
luminosity distance by 



Mz) = 51og 



diiz) 



1 Mpc 



+ 25 = Slog 



Dl{z) 



+ /^o , (12) 



where Di{z) = Hodiiz) jc is the Hubble free luminosity dis- 
tance and [iQ = 42.38-51ogjQ/i, with h being a dimensionless 
parameter defining the value of the Hubble parameter at the 
present epoch as Ho = 100/z km s"' Mpc"'. Using the para- 
metric form of the di (Eq. the distance modulus can be 
expressed in terms of the parameters a and b as 



^th(fl,^;z) = 5iog 



z(l+az) 
l+bz 



(13) 



The observed values of the distance modulus /iobs(Zi) corre- 
sponding to measured redshifts Zi are given in terms of the 
absolute magnitude M and the apparent magnitudes niobsizi) 
by 



l^ohs{Zi) = mobAzi)-M . 



(14) 



To obtain the best-fit values of the parameters a and b from 
SNe la observations we perform a analysis which involves 
minimization of suitably chosen function with respect to 
the parameters a and b. For our analysis of SNe la data we 
use the x^ fu nction consider ed in (Xu et al. 2010). We refer 
the reader to (IXu et al. 20101) for a comprehensive discussion 
on the choice of x^ function and its minimization. The x" 
function (for the analysis of SNe la data) is first defined in 
terms of parameters a, b andM' = fio+M (called the nuisance 
parameter) as 



(a,b,M') = Yl 



(^iobsiz^- l^ih{a,b,z)Y 



i=\ 

N 



^ (5 logio(D£(a, fe, z,)) - mpbsizi) + M'f^ ^^^ 



/=1 



where cr, is the uncertainty in observed distance mod- 
ulus and is the total number of data points. Its 
marginalization over the nuisance parameter as Xsn^'^'^) ~ 



-21n/_'^exp[-|x^(a,^,A/')]'iM' leads to ^^^(a,h) = P- 
(27/?) + ln(/?/27r) with 



(5 logio(Dz,(ci, ^, z;)) - niohsizdf 



(=1 

N 



g ^ ^ (5 logio(Dz.(a, /?,z,)) - mohs(zd) 



" 1 



/=i 



The function ~x^{a^b,M') has a minimum at M' = Q/R which 
gives the corresponding value of h as 10*'^"*' +42.38)/5 jjj-Qp. 
ping the constant term ln(/?/27r) from Xsw function 



xl^{a,b) = P-^ 



(16) 



can be used in for the likelihood analysis. 

The observation of baryon acoustic oscillations provides 
another evidence for the existence of dark energy. In the early 
universe, free electrons and protons were coupled with highly 
energetic photons of the relativistic plasma through scatter- 
ing. The high pressure in the plasma drives the primordial 
cosmological fluctuations to propagate outward at a relativis- 
tic speed. After the universe had cooled down sufficiently, at 
a certain point, electrons and protons combine to form neu- 
tral hydrogen atoms and thereby decoupling photons from 
baryons. This causes an abrupt fall of the speed of prop- 
agation of the acoustic wave. The baryon acoustic oscilla- 
tion leaves their signature on the large scale structure of the 
universe. The Slogan Digital Sky Survey (SDSS) measures 
the correlation function of the large sample of luminous red 
galaxies. The acoustic peak detected by them provides a stan- 
dard ruler with which the absolute distance corresponding to 
a typical redshift z = 0.35 can be determined. The standard 
ruler is given by the dimensionless parameter A in Eq. (|5]l as 
A = 0.469 ±0.017. The dimensionless quantity A is as given 
in Eq. (|5]l is constructed from the following set of equations 



A{z) = Dv{z) 



Dviz) -- 



Dl(z)z 



-I 1/3 



H{z) 



,Da(z) = 



dz' 



Hiz') 



^17) 



where Da(z) is the comoving angular diameter distance and 
Dviz) is the dilation scale. Considering a, b and as param- 
eters to be determined from observations, the function for 
analysis of BAO data is given as 



2 /oO 
ABAO'."m 



a,b) = 



[Ainl.aM-A,^,^ 
(AA)2 



(18) 



with Aobs =0.469, ZU = 0.017 and A(17° ,a,fo) being given by 
Eq. i). 

Determination of Hubble parameter from observational 
measurements is another probe to the accelerated expansion 
of the universe attributed to the dark energy. Compilation 
of the observational data based on measurement of differ- 
ential a ges of the galaxies b y Gemini Deep Deep Survey 
GDDS ^Abraham et al. 20041) . SPICES and VDSS surveys 
provide the values of the Hubble paramet er at 15 differ- 
ent redshift valu es (Gaztanaga ef a/. 2009; Ri ess et al. 20091: 
IStern et al. 20101; ISimon et al. 2 005). The x^ function for the 
analysis of this observational Hubble data can be defined as 



15 r 



Xohd(«.^) = X1 



H{a ,b;zi)-Hoi„{zi) 



(19) 
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where Hobs is the observed Hubble parameter value at z, with 
uncertainty E,-. 

Varying the parameters a, b and $7° freely we minimize the 
function which is defined as 

x\a, bMl) = xln{a,b) + Xbao(«, ^, ^l) + Xohd(«, ^) -(20) 

The values of the parameters a, b and Vl^^ at which minimum 
of is obtained are the best-fit values of these parameters for 
the combined analysis of the observational data from SNe la, 
BAO and OHD. With these values of the parameters we find 
the variation of the dark energy equation of state wx{z) and 
the dark energy density parameter (ilx(z) = 1 - ^miz)) using 
Eqs. (|9]l and (fTTT l respectively. We also find the la ranges of 
the parameters a, b and O", from the analysis of the observa- 
tional data discussed above. In this case of three parameter 
fit, the 1(7 (68% confidence level) allowed ranges of the pa- 
rameters correspond to < Xmin + ^X^^ where Ax^(= 3.53) 
denotes the la spread in x^ corresponding to three parame- 
ters. For these allowed domains of the parameters, we also 
obtain the corresponding la ranges of the quantities wx{z), 
Vlxiz) and r2,„(z) from Eqs. (|9]l and (fTTl i. The computation of 
Vt„,{z) from Eq. (fTTT l with the parameters a, b and O", as inputs 
from their 1 a ranges obtained in a way described above does 
not directly ensure that the condition Vl„,(z) < 1, which fol- 
lows from the definition of Vl„,{z), is always respected. (How- 
ever, we have seen that for the best-fit values of parameters a, 
b and 57° as obtained from the analysis, Vl„,(z) lies below 1 
for the range of z probed by SNe la observations.) To circum- 
vent this, VLmiz) at some particular value of z corresponding 
to an earlier epoch (beyond the range of measured redshifts 
z ^ 1 .76 of SNe la events) is not allowed to exceed some cho- 
sen benchmark value (say, a) below 1 . To take into account 
this constraint we find the domain of the (a,fe, 57„,) parameter 
space for which the two conditions, viz. x^ ^ Xmin + ™d 
r2,„(a,Z7,r25),;z) < a are simultaneously satisfied. The Icr range 
of the parameters thus obtained are, therefore, dependent on 
the initial condition of matter density at some earlier epoch 
which we choose here as z = 2. We study the impact of im- 
posing the constraint {Vl„,{z = 2) < a) on the parameter space 
by finding allowed domains of the parameter space for differ- 
ent choices of values of a. We also study this effect on the 
allowed region of the equation of state (yvxiz)) of dark energy. 

4. RESULTS AND DISCUSSIONS 

In this section we present and discuss the results of 
combined analysis of the SNe la, OHD and BAO data. As 
mentioned earlier we have considered five different sets of 
data for the analysis of observations of Supernova la along 
with the OHD and BAO data. The SNe l a data sets con- 
sidered here are HST-i- SNLS +ESSENC E (Riess et al. 20071: 
IWoodVasey et al. 20071: iDavis et al . 2007), SALT2 
data an d MLCS data (iKessler ef a/. 20 09). UNION 
data dKowalski et al 20081 and UN10N2 data 
(lAmanullah et al. 20101) . In the analysis, we have taken 
each one of these five sets of SNe la data at a time with 
OHD and BAO data to compute x^ = Xsn + Xohd + Xbao 
different sets of values of the parameters a, b and 51° . 

In Table [U we present the best-fit values of the parameters 
a, b and 57° obtained from analysis of different data sets. The 
minimum value of x^ along with total number of data points 
for each set of data are also shown. We refer to the 5 different 
sets of data considered here by 'set I', 'set IT, 'set III', 'set IV' 
and 'set V as accordingly listed in column 1 of Table [T] We 



SNe la data sets best-fit values of Minimum 



+ BAO + OHD (a, h, n° ) value of 



HST+SNLS+ESSENCE 
+ BAO + OHD 
(Data Set: I) 
No. of data points = 192+1+15 


(1.437,0.550, 0.268) 


199.267 


SALT2 
+ BAO + OHD 
(Data Set: II) 
No. of data points = 288+1+15 


(1.401, 0.542, 0.272) 


560.083 


MCLS 
+ BAO + OHD 
(Data Set: III) 
No. of data points = 288+1+15 


(1.401, 0.653, 0.296) 


783.078 


UNION 
+ BAO + OHD 
(Data Set: IV) 
No. of data points = 307+1+15 


(1.635,0.699, 0.268 ) 


311.615 


UNION2 
+ BAO + OHD 
(Data Set: V) 
No. of data points = 557+1+15 


(1.289, 0.458,0.272) 


544.074 



TABLE 1 

Best-fit values of parameters and minimum values of for 

EACH of the 5 DATA SETS CONSIDERED. 



then find the Icr ranges of the parameters a, b and 57° for each 
of the data sets (I-V) using the method described in Sec.|3] To 
ensure the boundedness: 57,„(z) < 1 , we obtain the allowed 
ranges of the parameters for different values of a (defined in 
Sec. O. The la allowed region of the parameters obtained 
from the combined analysis of SNe la, OHD and BAO data 
are presented in the planes of any two parameters of set { a, b, 
57° } by marginalizing over the third one. In Figures [T] |2]and 
[3] we show the la contours in parameter planes a-b, a -57° 
and 57° respectively for different choices of the values of a 
(57,„(z = 2) < a). The five panels (row-wise) from left to right 
in each of the Figures (Figure[T]|2]and|3]l correspond to results 
of analysis of data sets I to V. For panels 1-4 (corresponding to 
data sets I-IV) in these three figures, the la contours are plot- 
ted for three different values of a viz. 0.8, 0.9 and 1.0 i.e. for 
57,„(z = 2) < 0.8,0.9 and 1. In the last panel corresponding to 
the data set V (UNION2 data along with OHD and BAO), the 
same is plotted for a = 0.9,0.95 and 1. The values of a below 
0.9 are not chosen for data set V (UNION2+OHD+BAO) be- 
cause the value of 57,„(z = 2) exceeds 0.9 even when calculated 
at the best-fit values of parameters a, b and 57° obtained from 
the analysis of data set V. The UNION2 data (along with OHD 
and BAO) thus restricts the matter density parameter value at 
an epoch z = 2 to lie slightly below 0.95 (at la level). 

With the best-fit values of the parameters a, b and 57° as 
obtained above (listed in Table [T] for different data sets) we 
compute the equation of state w(z) of dark energy as a function 
of redshift z using Eq. The plots for w(z) vs z are shown 
by solid curves in Figure|4]for all the 5 data sets. Plots from 1 - 
5 (row-wise) correspond to data sets I-V. Using the 1 a range 
for the parameters a,b and 57° as obtained above for different 
values of a (flmiz) < a) we obtain the corresponding spread 
in w(z). The shaded regions in Figure|4]show these la bands 
of w(z). This has been shown for different choices of the value 
of a. The first 4 plots from left to right in the upper (lower) 
panel are for data sets I-IV respectively with 57,„(z = 2) < 0.8 
(57m(z = 2) < 0.9). The two plots in the extreme right column 
correspond to the data set V with the two constraints: 57m(z = 




Fig. 1. — 1(T contour in the a - h parameter space (with marginalization 
over nj,?') obtained from analysis of data sets I-V (see text) from left to right. 
In each of the 4 figures from left, contours are shown for three different values 
of a = 0.8, 0.9 and 1 . The figure in the extreme right (for data set V) contours 
are shown for a = 0.9,0.95 and 1 (described in text). 
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Fig. 2. — l(j contour in the a - parameter space (with marginahzation 
over b) obtained from analysis of data sets I-V (see text) from left to right for 
same set of values of a as in Figure ^ 
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Fig. 3. — 1(T contour in the b - Slj), parameter space (with marginalization 
over a) obtained from analysis of data sets I-V (see text) from left to right for 
same set of values of a as in Figure[T] 



Fig. 4. — Plots of wx(z) vs z for different data sets : Columns from left 
to right con'espond to data sets I-V respectively. Figures for data set I-IV are 
plotted for Qm(z = 2) < 0.8 (upper panel) and f2m(z = 2) < 0.9 (lower panel). 
For data set V (UNION2). plots are shown for n„(z = 2) < 0.9 and 0.95 at 
upper and lower panels respectively. 
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Fig. 5. — Plots of Qx(z) and Qm(z) vs z for different data sets : Columns 
from left to right correspond to data sets I-V respectively. Figures for data 
set I-IV are plotted for n,„(z = 2) < 0.8 (upper panel) and n„,(z = 2) < 0.9 
(lower panel). For data set V (UNION2), plots are shown for C„,(z = 2) < 0.9 
and 0.95 at upper and lower panels respectively. 



2) < 0.9 (upper panel) and Qm(z = 2) < 0.95 (lower panel). 
The same best-fit curves for wxiz) vs z are plotted both in the 
upper and lower panels of a given column. From Figure 2]we 
observe that, in some cases, the wxiz) vs z plots (solid curves) 
corresponding to the best-fit values of the parameters a, b and 
fijj, lie well within the respective Icr regions. They barely 
remain within such regions in some other cases. In one other 
occasion, the best-fit plot is outside the Icr region for most of 
the range of z considered. For example, for data sets III and IV 
the best-fit curve remain within the corresponding 1 a range 
for both the choices: ^^^(z = 2) < 0.8 (upper panels, columns 
3 & 4) and ^^(z = 2) < 0.9 (lower panel, columns 3 & 4). 
For data sets I and II as shown respectively in first and second 
columns of Figure |4|i, however, one can see that although the 
best-fit lines are within the Ict regions obtained for ilmiz = 
2) < 0.9 but certain segments of the best-fit lines tend to come 
out of the respective Icr domains obtained for the case ^miz = 
2) < 0.8. In the last column of Figure|4] corresponding to data 
set V, we show the Icr spread in the variation of wxiz) with z 
for the two choices: f2m(z) < 0.9 (upper panel) and f2m(z) < 
0.95 (lower panel). As evident from the plots, a large segment 
of the best-fit wx{z)-z curve in this case is not contained within 
the Icr region obtained for the case 51„,(z = 2) < 0.9 where as 
fully contained within the Icr region obtained for Vl,n{z = 2) < 
0.95. 

The above results can be interpreted and summarized like 



this: The data sets I ((HSTh-SNLSh-ESSENCE)h-BAOh-OHD) 
and II (SALT2h-BAOh-OHD) support the fact that matter den- 
sity parameter Q,„ at an early stage of the universe at z = 2 
was ^ 0.8 whereas the data sets III (MLCSh-BAOh-OHD) and 
IV (UNIONh-BAOh-OHD) can accommodate values of matter 
density parameter at the epoch z = 2 even a bit lower than 0.8. 
According to analysis of data set V (UNION2h-BAOh-OHD), 
the matter density parameter at z = 2 is only allowed to have 
values greater than 0.9. 

Another important observation in the context of Figure 
ID is the following. The variation of wxiz) over the range 
^ z ^ 1.76 (relevant for SNe la) with its Icr spread obtained 
from data sets I, II and IV with VL,„{z = 2) < 0.8 (columns 1, 
2 and 4 in the upper panel of Figure |4] show non-overlap be- 
tween the Icr bands of wxiz) at two well separated z values 
in the range ^ z ^ 1 .76. This phenomenon points to obser- 
vational evidence (at Icr level) for varying equation of state 
of dark energy as opposed to cosmological constant solution. 
On the contrary, such signatures of varying dark energy are 
not obtained from the analysis of the sets of data I, II, IV 
and V with Vl,„{z = 2) < 0.9. This is evident from the plots 
in the lower panels of columns 1, 2 and 4 and in the upper 
panel of column 5 of Figure|4] The other data set III also does 
not provide any signature of varying dark energy (at 1 cr level) 
even for 0,„(z = 2) < 0.8. Therefore, knowledge of the matter 
density of the universe at some earlier epoch is instrumental 
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in providing observational evidences in favor of varying dark 
energy or cosmological constant solutions. 

Using Eq. ( fTTT i we compute the matter density parameter 
flmiz) as a function of z for the best-fit values of the param- 
eters a, b and ilj],. The Icr range of the parameter ri,„(z) for 
different choices of the value of a are also computed. The 
results are shown in Figure |5] by solid lines and shaded re- 
gions respectively. The plots for dark energy density parame- 
ter Vlx(z) = 1 -ilm(z) are also shown. Columns 1 to 5 respec- 
tively correspond to data sets I-V. For the upper and lower 
panels of columns 1-4, imposed constraints on Vl,„{z = 2) are 
= 2) < 0.8 and Vl„,{z = 2) < 0.9 respectively. For col- 
umn 5 plots are shown for r2„,(z = 2) < 0.9 (upper panel) and 
17m(z = 2) < 0.95 respectively. The dark energy starts domi- 
nating at the value of z where plots for il„,(z) and Vlx{z) in- 
tersect. The estimations of such an epoch for each of the data 
sets considered may be obtained from the point of intersec- 
tions of the corresponding best-fit plots of Vl„,{z) and ^xiz)- 
For example, this intersection point are obtained at z ~ 0.4 
for data sets IV and V whereas for data sets III the 2 plots 
intersect at z — 0.48. As in Figure HI in this figure also, the 
best-fit plots for data set I and II for the choice ri,„(z = 2) = 0.8 
do not remain contained within the corresponding 1 a regions. 
The same is true for the data set V (upper panel of the last 
column of Fig. 5; rim(z = 2) = 0.9). From Figure |5] we also 
obtain the 1 cr uncertainty for ilJJ, - the matter density param- 
eter at present epoch. For data sets I, II, IV and V the Icr 
spread of iljj, is ~ 0.23-0.31 and for data set III this spread 
is ^ 0.25-0.35. These ranges are calculated with the choice 
r2„,(z = 2) = 0.8 for data sets I,II,III and IV while for data set 
V the estimation is for the choice J7m(z = 2) = 0.95. We like 
to comment that these ranges will be affected negligibly for 
other choices of ri,„(z = 2) considered here as is obvious from 
Figures. |2]and[3] The earlier discussed facts that the data sets I 
and II allow matter density parameter values at an early stage 
of the universe at z = 2 <; 0.8 whereas data sets III and IV 
can accommodate even lower values of matter density param- 
eter at the epoch z = 2 and data set V only allows values a 
bit greater than 0.9 are also reflected in the plots of Figure 
|5] through the containment or non-containment of the best-fit 
plots within the corresponding 1 tr limits. 

5. CONCLUSION 

In this paper we have performed a combined analysis 
of the SNe la, OHD and BAO data assuming a closed 



form parametrization of the luminosity distance in terms 
of redshift. Five different sets of SNe la data designated 
as (HSTh-SNLSh-ESSENCE), SALT2, MLCS, UNION and 
UNION2 are independently analysed combining the individ- 
ual data sets with BAO and OHD. Among the 5 SNe la data 
sets, the UNION2 data is the most recent. From the analysis 
we find the best-fit values of the parameters a and b appear- 
ing in the luminosity distance - redshift parametrization along 
with the matter density parameter at the present epoch (ilj],). 
For this we use a 3-parameter fit to each of the SNe la data 
sets along with OHD and BAO. We also find the Icr ranges 
of the parameters a, b and imposing a constraint on the 
matter density parameter (fi,,,) that at earlier epoch (chosen 
here as z = 2) its value does not exceed a certain value a. Re- 
sults are presented for some benchmark values of a as 0.8, 
0.9, 0.95 etc. We also compute z-variations of the equation of 
state of dark energy {wxiz)) and matter and dark energy den- 
sity parameters {Cl,r,{z) and ^xiz)) corresponding to the best- 
fit values of the parameters a, b and iljj, and their Icr ranges. 
The results of the analysis show that knowledge of the matter 
density of the universe at some earlier epoch is instrumental 
in providing observational evidences in favor of varying dark 
energy or cosmological constant solutions. Also the SNe la 
data sets (HSTh-SNLSh-ESSENCE) and SALT2 (along with 
BAO and OHD) restrict the matter density parameter value 
at the earlier epoch (z = 2) not to go below 0.8. The data 
sets MLCS and UNION 1, on the other hand restrict the same, 
above a value which is slightly greater than 0.8. The recent 
UNION2 data, however, constrains the value to lie always 
above 0.9. We obtain the matter density parameter at present 
epoch, 9P„^ - 0.27 for four sets of data (I, II, III and V) while 
for data set IV, ri,°„ ~ 0.3. The Icr spread for ilj], lies within 
the range ~ 0.23-0.31 for data sets I, II, IV and V whereas 
for the data set III this range is ^ 0.25-0.35. The nature of 
variations of the dark energy equation of state are similar for 
data sets I-IV and is different for the data set V. The analysis 
of the data sets I, II, IV and V show that the dark energy starts 
dominating the matter from the epoch z ^ 0.4 and the same 
from the analysis of data set V is found to be z ^ 0.48. 



ACKNOWLEDGEMENTS 

The authors thank Anjan Ananda Sen for some useful dis- 
cussions. 



REFERENCES 



Abraham R. G. et al. , 2004, Astron. J. 127 2455 
Amanullah R. et al, 2010, Astrophys. J. 716, 712 
Astier R, 2000 arXiv : astro-ph/000830 6 
Chiba T.& Nakamura T., 2000, Phys. Rev. D 62 121301 
Chevallier M. & Polarski D., 2001, Int. J. Mod. Phy. D. 10, 213 
Corasaniti R S. & Copeland E. J., 2003 Phys. Rev. D 67 063521 
Daly R. A. & Djorgovsky S. G., 2003, Astroph. J. 597 9 
Davis T. M.et al., 2007, Astrophys. J. 666, 716 
Eisenstein D. J. et al. , 2005, Astrophys. J. 633, 560 
Fay S. & Tavakol R. , 2006, Phys. Rev. D 74 083513 

Gaztanaga E., Cabre A. & Hui L., 2009, Mon. Not. Roy. Astron. Soc. 399, 
1663 

Gerke B. & Efstathiou G., 2002 Mon. Not. Roy. Ast. Soc. 335 33 
Gong Y., 2005 Int. J. Mod. Phys. D 14 599 

Guo Z. K., Ohta N. & Zhang Y. Z., 2005 Phys. Rev. D 72 023504 

Huterer D. & Turner M.S., 1999, Phys. Rev. D 60 081301 

Huterer D. & Cooray A. , 2005, Phys. Rev. D 71 023506 

Jassal H. K., Bagla J. S. & Padmanabhan T., 2005 Mon. Not. Roy. Astron. 

Soc. 356 LI 1 
Jimenez R. & Loeb A., 2002, Astrophys. J. 573, 37 



KesslerR. etal, 2009, Astrophys. J. Suppl. 185 , 32 

Komatsu E. et al. [WMAP Collaboration], 201 1, Astrophys. J. Suppl. 192, 18 
Kowalski M. et al., 2008, Astrophys. J. 686, 749 
LinderE. V., 2003 Phys. Rev. Lett. 90 091301 
Ma C. & Zhang T-J., 201 1, Astrophys. J. 730, 74 

Maor I., Brustein R., McMahon J. & Steinhardt R J., 2002, Phys. Rev. D 65 
123003 

Nesseris S. & Perivolaroupolos L., 2004 Phys. Rev. D 70 043531 
Perlmutter S. et al. , 1999, Astrophys. J. 517, 565 
Riess A. G. et al. , 1998, Astron. J. 116, 1009 
Riess A. G.etai, 2007, Astrophys. J. 659, 98 
Riess A. G.etai., 2009, Astrophys. J. 699, 539 

Roy Choudhury T. & Padmanabhan T., 2005, Astron. Astrophys. 429 807 
Saini T. D. , 2003, Mon. Not. Roy. Ast. Soc. 344 129 

Sahni V., Saini T. D., Starobinsky A. A. & Alam U., 2003 JETP lett. 77 201 
Shafieloo A. , 2007, Mon. Not. R. Astron. Soc. 380 1573 
Simon J., Verde L. & Jimenez R., 2005, Phys. Rev. D 71 123001 
Starobinsky A, 1998, JETP Lett. 68, 757 

Stem D., Jimenez R., Verde L., Kamionkowski M. & Stanford S. A., 2010, 
JCAP 1002, 008 



Wang Y.& Lovelace G., 2001, Astroph. J. 562 LI 15 
Wang Y. & Mukherjee P., 2004, Astroph. J. 606 654 
Weller J. & Albrecht A., 2002, Phys. Rev. D 65, 103512 
Wetterich C, 2004, Phys. Lett. B 594 17 
Wood-Vasey W. M. et al, 2007, Astrophys. J. 666, 694 



7 



Wu P & Yu H., 2006, Phys. Lett. B 643 315 
Xu L. & Wang Y., 2010 JCAP 1006, 002 
Zhang Q-J. & Wu Y-L., 2010, JCAP 1008, 038 



